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Resumen 


Se resuelve de manera exacta la ecuación Fokker-Planck con un término de 
sumidero, utilizando la difusividad de Fujita y la relación de Parlange entre la 
conductividad y difusividad. Para obtener la solución, se introduce primero 
el potencial de Kirchhoff y enseguida la coordenada móvil de Fujita-Storm, la 
ecuación diferencial toma la forma de la ecuación de Burgers, que es lineal en 
el término difusivo. El coeficiente convectivo de esta última se sustituye por 
la transformación de Hopf-Cole, con la finalidad de deducir la ecuación lineal 
clásica del calor. Durante las transformaciones, el término de sumidero se define 
funcionalmente, de modo que el resultado final sea precisamente la ecuación de 
calor sin término de sumidero. La solución exacta del potencial de Hopf-Cole 
se obtiene con la transformada clásica de Laplace para algunas condiciones 
iniciales y de frontera de interés. La solución de la ecuación Fokker-Planck 
en el espacio físico se obtiene a través de la inversión de las transformaciones 
utilizadas. La solución incluye como casos particulares las soluciones de Sanders 
et al., y de Broadbridge y White. La solución exacta puede ser utilizada para 
validar soluciones numéricas de la ecuación Fokker-Planck y en estudios sobre la 
extracción de agua por las raíces de las plantas. 


Palabras clave: características de Fujita-Parlange, transformación de Kirchhoff, 
transformación de Hopf-Cole, ecuación de Burgers. 


Introducción y de la ley de Darcy generalizada a los suelos no 
saturados: 
Considérese la ecuación Fokker-Planck del 00 
movimiento unidimensional del agua en el q=-D(0) ed K(0) (3) 
suelo: 
en donde 0 es el contenido volumétrico de agua, 
E = , 0 E -Y(2,t) (1) también denominado contenido de humedad; 
dt oz dz | d8 dz 


q, el caudal de agua por unidad de superficie 
de suelo o flujo de Darcy; z, la coordenada 


que resulta de la combinación de la ecuación 
de continuidad implicada por el principio de 
conservación de la masa: 


Ylz,t) (2) 


vertical orientada positivamente hacia abajo; t, 
el tiempo; Y, un término de sumidero, como el 
volumen de agua extraído por las plantas por 
unidad de volumen de suelo en la unidad de 
tiempo; K(8), la conductividad hidráulica no 
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saturada, y D(0) es la difusividad hidráulica 
definida por: 


D(O)K 0) 4) 


en la cual y(0) es la presión del agua en el 
suelo como una función del contenido de 
humedad o curva de retención de humedad, 
y es expresada como la altura de una 
columna equivalente de agua. K(0), D(0) y 
y(8) son conocidas como las características 
hidrodinámicas del suelo. 

Se han construido pocas soluciones ana- 
líticas exactas de la ecuación Fokker-Planck 
para características hidrodinámicas no lineales. 
Entre ellas está la solución en serie en la raíz 
cuadrada del tiempo presentada por Philip 
(1957) para características hidrodinámicas 
más o menos arbitrarias, y condiciones límites 
específicas, válida para tiempos cortos y de 
radio de convergencia desconocido. 

También se cuenta con soluciones analí- 
ticas para condiciones límites específicas, 
válidas para todo tiempo, y que sonexactas para 
características hidrodinámicas particulares, 
como la solución lineal (e.g. Philip, 1969); la 
solución de Knigth (1973), adaptada de una de 
las soluciones de la ecuación de Burgers (1948); 
la solución cuasi-lineal de Fuentes et al. (2001), 
que contiene como casos particulares las dos 
soluciones precedentes. En la solución cuasi- 
lineal se asume una difusividad hidráulica 
constante y una conductividad hidráulica, que 
resulta de una interpolación entre la variación 
lineal de la conductividad con respecto del 
contenido de humedad (suelo lineal) y de 
la variación cuadrática en relación con este 
contenido de humedad (suelo de Knigth). 

Por otro lado, se cuenta con la solución 
analítica de Parlange et al. (1982), que es 
exacta cuando la difusividad es asimilada a 
una densidad de Dirac, y la conductividad 
representa, por un lado, una función lineal con 
respecto al contenido de humedad y, por el otro, 
cuando la derivada de la conductividad con 
respecto al contenido de humedad es también 


asimilada a una densidad de Dirac. El primer 
caso corresponde a la solución clásica de Green 
y Ampt (1911), y el segundo a la solución de 
Talsma y Parlange (1972). 

Fujita (1952) ha deducido una solución 
no lineal de la infiltración unidimensional 
horizontal, es decir, del movimiento del agua 
en ausencia de gradientes gravitacionales, 
conocido como absorción, representada por 
la ecuación (1), sin los términos segundo 
y tercero del miembro derecho 6, = [D(8)8_], 
donde los subíndices indican las variables 
con respecto de las cuales se lleva a cabo 
la diferenciación; la difusividad no lineal 
de Fujita contiene como casos particulares 
los comportamientos extremos, es decir, la 
difusividad constante y la difusividad de 
Dirac. 

Bluman y Kumei (1980) han demostrado 
que la ecuación de la absorción es invariante 
bajo transformaciones de Lie-Bácklund si y 
sólo si la difusividad es precisamente la de 
Fujita, y que cualquier ecuación diferencial 
invariante bajo estas transformaciones puede 
ser asociada con la correspondiente ecuación 
diferencial parcial lineal. Fokas y Yortsos (1982) 
han extendido estos resultados a la ecuación de 
flujo bifásico de la absorción 6, = [D(8)8_] — v(t) 
f.(0), y deducen la solución para una condición 
de flujo constante; la función de fracción de 
flujo f(0) resulta proporcional a la integral 
de la difusividad de Fujita y v(t), el flujo 
combinado de ambos fluidos considerados 
incompresibles, es tomado constante. Esta 
solución formalmente corresponde a la solución 
para flujo constante de la ecuación (1), con Y = 
0, cuando dK/d8 « D(0), es decir, cuando K(y) 
es una función exponencial (Gardner, 1958). 

En esta misma dirección, Rogers et al. 
(1983) extienden los resultados para deducir 
la solución para flujo constante de la ecuación 
de flujo bifásico en un campo gravitacional 6, 
= [D(0)6_] — of.(0) — (K(8)[1-f(0)]).; los autores 
han asumido K(8) «< 8 y v constante. Sanders et 
al. (1988) han adaptado la solución de Rogers 
et al. (1983), como una solución de la ecuación 
(1) con Y =0, haciendo la equivalencia K'(0) 
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- of'(0) + (K(0)[1 — A(0))),, La nueva función 
de conductividad permite llevar la ecuación 
Fokker-Planck a la forma de la ecuación de 
Burgers (1948), la cual puede ser linealizada 
con una transformación de Bácklund; la 
función de conductividad obtenida, al igual 
que la difusividad de Fujita, es no lineal. 
Broadbridge y White (1988) también 
adaptan la solución de Rogers et al. (1983) y 
presentan un análisis de las características 
Esta sido 
extendida a condiciones iniciales de cierto 


hidrodinámicas. solución ha 
interés práctico y al drenaje de una columna 
de suelo (Warrick et al., 1990; Warrick y Parkin, 
1995; Parkin et al., 1995). La solución también 
se ha usado para validar soluciones cuasi- 
analíticas del fenómeno de la infiltración (Si y 
Kachanoski, 2000). 

Sin embargo, con las características hidro- 
dinámicas presentadas por Broadbridge y 
White (1988), no es posible representar, por 
ejemplo, el importante caso teórico extremo 
de Green y Ampt (1911); tampoco estas 
características contienen el suelo de Gardner 
(1958), de gran importancia teórica en el 
estudio del flujo monofásico multidimensional 
del agua en el suelo (Wooding, 1968). Estas 
características también inducen que la 
capacidad máxima de almacenamiento en el 
fenómeno de ascenso capilar no está acotada 
para todo tipo de medio. 

El presente trabajo tiene como objetivos: 


+ Demostrar que las características hidrodi- 
námicas resultantes del procedimiento 
de reducir la ecuación Fokker-Planck a la 
ecuación de Burgers corresponden a una 
relación entre la difusividad y la conduc- 
tividad propuesta por Parlange et al. (1982), 
si en ésta se introduce la difusividad de 
Fujita (1952). Las características de Fujita y 
Parlange contendrán como casos particu- 
lares las utilizadas en el establecimiento de 
las soluciones cuasi-lineal y generalizada 
de Dirac de la infiltración. 

e Establecer la forma requerida del término 
de sumidero en la ecuación Fokker-Planck, 


de modo que ésta se reduzca siempre a la 
ecuación de Burgers. 

e  Deducir de manera concisa la solución 
exacta de la ecuación Fokker-Planck, 
con las características hidrodinámicas de 
Fujita y Parlange, correspondiente a una 
condición de flujo constante en la superficie 
de una columna semi-infinita de suelo y a 
algunas condiciones de humedad inicial 
de interés. 


Las condiciones límite de la ecuación 
Fokker-Planck 


Para resolver la ecuación (1) es necesario 
proporcionar las condiciones inicial y de 
frontera, y las características hidrodinámicas. 
La condición inicial en una columna semi- 
infinita de suelo tendrá la forma: 
£=0, .0<z, 6=0(z) (5) 
En la frontera superior puede especificarse 
el contenido de humedad o el flujo de Darcy, 
a saber: 


(6) 
-D(0)7+K(0)= 8. 


donde 6... y q, pueden ser funciones del 
sup Ss 
tiempo. 
En la frontera inferior de una columna 
semi-infinita la condición es: 


t>0, 0=lim8,(z) (7) 


Las características hidrodinámicas se esta- 
blecen en la siguiente sección. 


Las características hidrodinámicas de 
Fujita y Parlange 


La difusividad hidráulica de Fujita (1952) 
puede ser escrita como (Fuentes et al., 1992): 
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K.-K - 
pre El le 1 0 5 A 
9,-0, (1-00) 
(8) 
6 DE 
0,—0, 


donde 6, es el contenido de humedad a 
saturación y K_ = K(0) es la conductividad a 
saturación; 9 <8 esun contenido de humedad 
de referencia y K, = K(0,) es la conductividad 
correspondiente; O es un grado efectivo de 
saturación; 0. es un parámetro adimensional; 
A. es la escala de Bouwer (1964), definida 
por: 


RE [D(oJa 


9 
j KK a, : 


La relación esencial entre la conductividad 
y la difusividad hidráulicas propuesta por 
Parlange et al. (1982) es la siguiente (Fuentes et 
al., 1992): 


K(0)-K A: 
=—=0|1-pB+pB/|D(0)40/ | D(O)48| (10 
OE -oli-0000(005/ [oy a 
donde f es un parámetro adimensional. 

La introducción de la ecuación (8) en la 


ecuación (10) proporciona la siguiente función 
de la conductividad K(O): 


o[1-P$+(P-a)0] 


(11) 


La función w(0) se obtiene mediante la 
sustitución de las ecuaciones (8) y (11) en la 
ecuación (4), y realizando la integración con 
la condición y(0 ) = y, una presión de burbuja 
en saturación que se asumirá, sin pérdida de 
generalidad, igual a cero. Cuando se toma 
9, = 0, donde 6, es el contenido de humedad 
residual definido, de modo que K(8) =0 y y(0) 
> —« (Brooks y Corey, 1964), se obtiene la curva 
de retención siguiente: 


vv | E. 


(1-a)0 
B 1-P+(P-0JO 0-0 ii 
—O0 —=B+(pB-0 _U—0, 
i B(1-B) .l (1-aJe , id O, Sl O, 
donde y, =-A.. 


Si se utiliza la notación sol(a, PB) para 
representar los diferentes suelos, el suelo 
cuasi-lineal corresponde a sol(0, B), es decir a 
una difusividad constante D = KA,/(8, -— 0) y 
a una conductividad de la forma K(O) = K.O[1 
- B + BO]; el suelo lineal es dado por sol(0, 0) y 
el suelo de Knigth (1973) por sol(0, 1). 

Cuando Aa > 1, la difusividad de Fujita 
se vuelve una densidad de Dirac. En efecto, 
haciendo O = 1 en la ecuación (8), la difusividad 
se concentra en saturación cuando oa. > 1; se 
tendrá D(0) = K 1,8(0 _—0). El suelo generalizado 
de Dirac es representado por sol(1,B); de la 
ecuación (10) se deduce para esta difusividad 
que dK/d0 = K [1 + BS(1 - 0)]; el suelo de Green 
y Ampt (1911) corresponde a sol(1, 0), y el suelo 
de Talsma y Parlange a sol(1, 1). 

El suelo de Broadbridge y White (1988) es 
representado por sol(a,, 1); la difusividad es la 
de Fujita y la conductividad es K(0) = K (1 — a.) 
0?/(1 - 010). La curva de retención se deduce 
de la ecuación (12): 


y (0) = y, fol ja 2 (13) 


Es claro que este suelo no puede represen- 
tar el suelo lineal ni el suelo de Green y Ampt, 
en los cuales $ =0. La solución no es tan versátil 


como han señalado los autores. 

El suelo de Gardner (1958) es representado 
por sol(a,, a), es decir, por la difusividad de 
Fujita y por la conductividad K(0) = K(1 
- 0)0/(1 - 940). Con la curva de retención 
correspondiente y(0) = wylin[(l - a0)/ 
(1 — 0JO], se deduce la función propuesta 
por Gardner K(w) = Kexpl-y/vy.), de 
importancia teórica en el estudio del flujo 
multidimensional del agua (Wooding, 1968) 
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y también de importancia práctica, ya que la 
solución de Wooding se usa para identificar 
parámetros del suelo a partir de pruebas de 
infiltración (Reynolds y Elrick, 1985). 

En el fenómeno del ascenso capilar 
de agua en una columna semi-infinita de 
suelo con un contenido de humedad inicial 
constante (0), las características de Fujita 
y Parlange proporcionan una capacidad de 
almacenamiento de agua en tiempos muy 
grandes acotada. La máxima capacidad de 
almacenamiento (M_) es dada por (Kunze et 


al., 1985): 


(14) 


Ms 9,)D(9) ¿y 


K(9)-K, 


La introducción de las ecuaciones (8) y (11) 
en la ecuación (14) proporciona la máxima 


capacidad de almacenamiento siguiente 
(Kunze et al., 1985): 
A.(0, —8 
M.= (0, mn ! (15) 
B 1-B 


Esclaroquelascaracterísticas de Broadbrigde 
y White proporcionan una capacidad de 
almacenamiento no acotada para todos los 
suelos ($ = 1), a diferencia de las características 
de Fujita y Parlange. En particular, para el 
suelo de Green y Ampt ($ = 0), la ecuación (15) 
proporciona M_ =1/(8,—0.). 

Así, 


que las características de Fujita y Parlange 


con este análisis se puede inferir 


son más versátiles que las de Broadbrigde y 
White. La pregunta que se plantea es que si 
las técnicas matemáticas —utilizadas para 
construir soluciones analíticas de la ecuación 
Fokker-Planck con las segundas— pueden ser 
aplicables en la construcción de soluciones 
con las primeras. En el apartado siguiente se 
demuestra que éste es el caso. 


Resolución de la ecuación Fokker-Planck 


La resolución analítica de la ecuación Fokker- 
Planck será grandemente simplificada con la 
introducción de las variables adimensionales, 
definidas como sigue: 


0-8  p__9-8 


9, = , ¡=> do  p, 
9, 0 (K,-K,)2. 
(16) 
ds K-K, 
K,-K, 
z K,-K, A 
E y Yi £ Y 17 
E RETO E 


Las ecuaciones (1), (2) y (3) permanecen 
invariantes; las variables tendrán asteriscos. 
Las ecuaciones (8) y (11) se transforman en: 


(18) 


D,. (9, 5] os 


9.[1-B+(P-0)0, ] 


K.(0,)= 1-08, 


(19) 


Para facilitar la escritura de las ecuaciones, 
se prescindirá de los asteriscos, en el entendido 
de que se trata de variables adimensionales 
hasta que no se diga lo contrario. 

Se introduce el potencial de Kirchhoff: 


(20) 


1 (0)= jo(a)aa 


el cual satisface u(1) = 1, para escribir la ecua- 
ción (3) como sigue: 


OS 1 + D(w)Y(z, A QD 


ou u 
=D 


La coordenada de Storm (1951) y Fujita 
(1952) definida por: 


(22) 


¿(21)=[[0(3.1)]az 


odian 
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permite linealizar el término difusivo de la 
ecuación (21), siempre y cuando f(0)0<1 ¿JD(0). 
Se demostrará que esta transformación 
linealiza la ecuación diferencial sólo en el caso 
de la difusividad de Fujita. 

De la ecuación (25) se deducen las reglas: 


6,6): 6) 


ES 


La ecuación (21) deviene: 


(23) 


(24) 


proa ot 


orefajea ala al Fle(z yo +Y(z, pez 
(25) 


en donde se ha utilizado la ecuación (2) y f'(0) 
= df/d6. El flujo de Darcy (q) está definido por 
la ecuación (3). 

Para linealizar el término difusivo se debe 
tener D(1)P(0) = k?, donde k es una constante. 
En virtud de las ecuaciones (2), (20) y (24), q 

-f(0)01/0€ + K(u); en consecuencia, para 
eliminar el cuadrado de la derivada en la 
ecuación (25), se debe asumir que df/d0 = -c,, 
donde c, es una constante. Se deduce f(0) =c,(1 
- 0.6), donde c, es una constante de integración 
y a. =C,/c, y por lo tanto D(0) = x?/c3(1 — 010). 
El potencial de Kirchhoff se deduce de la ecua- 
ción (20): (0) = x0/c2(1 — 00), y puesto que 
u(1) = 1, se tiene k?/c? = 1 - a. De este modo 
se obtiene D(8) = (1 — 0)/(1 — 010), la cual 
corresponde a la difusión de Fujita (1952), en 
escritura adimensional definida por la ecuación 
(18). Además: 


(1-a)0 


u(0)== 


o(u)=%— (26) 


1-0. +au 


r(0)- ol 
sa] 


1-0+au 


(27) 


169 =+| 


Con la ecuación (27), la ecuación (25) se 
simplifica a: 


gl) K d 
da a 0 o.+041)K (u)] 


97 yoon, RD) 


(28) 


a a 


a jr sa 


=K 


donde q,(t) = q(z =0,t) es el flujo de Darcy en la 
frontera superior. 

La ecuación (28) puede ser llevada a la 
forma de la ecuación de Burgers (1948) si los 
coeficientes del término convectivo tienen un 
comportamiento lineal en 4. Haciendo d[(1 - 
a. + 041)K(u)]/du = PB, + 28u, donde PB, y B son 
parámetros, se deduce: 


u(1-P+Pu) 


K = 
(4) 1-0 +au 


(29) 


en donde se han utilizado las condiciones K(0) 
=0 y K(1) = 1, que resultan en f, =1 -— P. 

La introducción de la ecuación (26) en la 
ecuación (29) conduce a la conductividad 
K(0) = 0[1 - $ + (PB — 098] /(1 — 010), la cual es 
precisamente la ecuación (19), que resulta de 
la relación entre la conductividad y difusivi- 
dad, propuesta por Parlange et al. (1982), 
cuando se introduce la difusividad de Fujita. 

Siguiendo con el mismo razonamiento, si 
la integral en la ecuación (28) tiene la forma 
o. — yu(z, t), entonces se debe tener o. = yu(0, 
t) para que la integral se anule en z = 0; y es 
un parámetro 


adimensional. Denotando 
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u (+) = (0, £), a la evolución del potencial de 
Kirchhoff en la superficie de la columna, se 
puede establecer: 


[rene =1[u(2.0)-u, (8)]. 


0 
Y (2,t)= a 


(0) 


es decir, en la ecuación (28) se tiene también: 


D(wY(Tt)= AAN a o)a (31) 


La ecuación (28) considerando las 
ecuaciones (29), (30) y (31), se transforma en la 


ecuación de Burgers (1948) siguiente: 


du _ en _ 


2 b— 32 
donde: 
_B-ya i-B-y(=0) 
E Yi : 
QL QL (83) 
0 
c= 
1-a 
o = 957 YM, (34) 


El parámetro k, siendo arbitrario, puede ser 
tomado como la unidad (Fuentes, 1992) o ser 
seleccionado de modo que la lámina de agua 
almacenada en la columna de suelo /(t) sea 
igual a la lámina almacenada en la “columna 
de Kirchhoff-Fujita”: 


(35) 


1(9)=j0(2,Nd2=fu(2 at 


en tal caso «=1/41-a. 

La ecuación de Burgers puede ser llevada 
a la forma de la ecuación lineal del calor, 
efectuando una sustitución no lineal debida 


a Hopf (1950) y Cole (1951), inicialmente 
propuesta para el caso de q, constante o c = 
0, restricciones no consideradas aquí. Esta 
transformación reemplaza el coeficiente del 
término convectivo de la ecuación (32), con 
la ayuda de una nueva variable (u), llamada 
potencial de Hopf-Cole: 


k(Qau+b-cq,)=-2k (36) 


2 9ln(u) 
ye 


esta transformación es un caso particular de 
las transformaciones de Bácklund. 
La ecuación (32) se transforma en: 


e 9 |1( du ¡2 uu _cdg, (87) 
05 | ul ot o? 2 dt 
de donde: 
du Yu c dq, 
lei [E+F(pJu (38) 


en la cual F(t) es una función desconocida del 
tiempo. 

El paso de la solución de la ecuación (38), 
u(€,t), a la solución de la ecuación (32), u(C,t), se 
realiza a través de la ecuación (36), a saber: 


=> je b-2x confu(E)]| (39) 


Para contar con la solución u(z,t) es nece- 
saria la relación entre z y €. Ésta se encuentra 
mediante la inversión de la ecuación (22), 
considerando las ecuaciones (27) y (39): 


6 
69-| pEnjz 
e (09, —b) 
2a (5 dE 


Cuando q, es una constante, la ecuación (38) 


Ll vas “(o 2) UN 


u(E1) 


se simplifica considerablemente, a saber: 
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2 
==> (41) 


la cual corresponde a la ecuación lineal del 
calor. 

De acuerdo con la ecuación (34), el flujo de 
Darcy en la frontera está dado por q,=q, + yu, en 
donde se ha supuesto que q, sea constante para 
deducir la ecuación (41). Si no hay término de 
sumidero (y = 0), las soluciones de la ecuación 
(41) serán aplicables para un flujo constante en 
la superficie (q, = q,), ya sea positivo (lluvia) 
o negativo (evaporación). Cuando existe el 
sumidero descrito por la ecuación (30) (y + 0), 
estas soluciones son aplicables al caso particular 
en que el flujo de Darcy en la superficie varíe 
como q,= 9, + YM. 

Para resolver la ecuación (41) es necesario 
transformar las condiciones límite. La con- 
dición de frontera en z = O es definida por 
la ecuaciones (6), (20) y (27): q,./D(u)= 


-x0u/08+K(u),/D(u). Obteniendo du/oí de 


la ecuación (36), utilizando las características 


hidrodinámicas, la propia ecuación (41) y 
considerando que uu (t) = u(0,f), se deduce: 


pt (22 1=0.- bo)a, 


+ (cg, ) di =0 


(42) 


La condición inicial se deduce de la propia 
ecuación (36): 


e a aa (E) +b- cq, ]u= 0 (43) 


Así, el problema consiste en resolver la 
ecuación (41), sujeta a la condición de frontera: 


u(0,t)=exp(Wt) (44) 
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y a la condición inicial: 


u(5,0)=exp [a (9)] (45) 
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ya que estas dos últimas ecuaciones satisfacen, 
respectivamente, las ecuaciones (42) y (43) con: 


dE y b2+2 (a Yl- a. bc) q, + (cq,) (46) 


A(C)== ant [loo 6)+)] (47) 


La resolución se vuelve simple, se 
descompone la función u en la suma de dos 
funciones que satisfacen la propia ecuación 
(41) (Carslaw y Jaeger, 1959): 


U=U +0, (48) 


La primera está sujeta a las condiciones 
inicial y de frontera: 


U (£,0)=0 ; U, (0,t)=exp (2+) (49) 
y la segunda está sujeta a las condiciones: 


u,(80)=exp[ A(8)]; u.(0,£)=0 (50) 


La solución del primer sistema se obtiene 
utilizando la transformada de Laplace y es 
(Carslaw y Jaeger, 1959, p. 64): 


O ad 


E : (51) 
O ido 7 PRIMICIA At 
a y 2xV/t 


donde erfc(x) = 1 — erfíx) es la función comple- 
mentaria de los errores y erf(x) es la función de 
los errores (Abramowitz y Stegun, 1970). 

La solución del segundo sistema (Carslaw y 
Jaeger, 1959, p. 59) se expresa como: 


E al 


u, (5 a) exp|— 


| at (52) 
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De las ecuaciones (51) y (52) se obtiene: 


du Do ant É 

E MAA o | 3 

oí 2k dl 2x 4/4 
A Ocenr E 

¿ep erfcl =—-A lt (53) 
2K J 2x4 


A ex E 
rada INTRO 
du 


C 1 í E 7 
0% = ao (E el 


ass - (E sel e 


4x4 


Eos 


4x4 


(54) 


Considerando las ecuaciones (26), (39), (40), 
(44), (48), (51), (52), (53) y (54), la solución está 
formada por: 


u(5,t) 


1-0 +0p(E,t) 2 


0(5,t)= 


160 Ze. 


2 pueo aca] E 
E] E de 


200) a Dl 


(57) 
ón exp (+) 
a [4 (E,t)+u, (Et) 


La solución de la ecuación Fokker-Planck 0 
= 0(z,t) tiene la forma paramétrica O = O(€,t) y 
z = z((,t), donde ( es el parámetro y q, el flujo 
efectivo. 

Una solución cerrada puede ser obtenida 
cuando A(í) presenta un comportamiento 


lineal, para ello es suficiente considerar que 
la función 0/(z) en la ecuación (5) es una 
constante. 

La solución incluye el caso q, > 0 (lluvia), 
9/2) una constante y $ = 1, presentado por 
Broadbridge y White (1988). Para la condición 
inicial, estudiada para PB = 1, propuesta por 
Warrick et al. (1990): 


(58) 


permite integrar la ecuación (52) de manera 
cerrada. 

La coordenada de Fujita €,, correspondiente 
a z,, se deduce de la ecuación (22): 


1-08, 
y1-oL 


E¡=Xx Z (59) 


La función A(€), considerando la ecuación 
(47), queda definida por: 


_ (41 /x)8 ¡O0<E€<E, 
e 3 pd 


donde: 


1 
Aj => [e - (20 pu, +b)); 


(61) 
A, => [cq, (204, +b)] 


en las cuales los potenciales de Kirchhoff Lu, 
y Hu, corresponden, respectivamente, a 0, y 
9, y son calculados con la ecuación (26). La 
función A(S) incluye el caso de un contenido 
de humedad constante en toda la columna (8); 
es suficiente hacer z, = 0 y 0, = 6. 

La introducción de la ecuación (60) en la 
ecuación (52) permite encontrar la función 
u (Ct) siguiente: 
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(ca) == exp[(a/)6+ 4%] pre: Ss +4) 


«/£ 


e 57 +a06) 
ep[-(A/06+ 4] e da Ayvi ) 
come ¿Eq pad 


7 azi] 


+(42/x)É + Ajtherfo ES 


xp J[(4,- 42, )/x]6 (42/10) 


+ Añtperfo ES Az yt | (62) 


y de aquí se deduce la derivada du. 108: 


STA Jer e j 


Ñ ol +40) 


Aj 2 6+61 
exp[- (4,/x)5 + Az der (ES 


lira] 


_ 5) 


op [ (A, AD/ JE + (2/0) (63) 
+ Añtherfe | a As Y 
Tex (A, A)/x JE (42/0% 


ambas expresiones necesarias para contar con 
la solución exacta planteada en las ecuaciones 
(55), (56) y (57). 

La función de sumidero se obtiene de la 
ecuación (34), considerando las ecuaciones (26) 
y (31), a saber: 


Via ot) 


(64) 
1-0a+opu(E,t) 0 


Y(E,+)=-yx 


« vj30]0u32], 


La derivada se obtiene de la ecuación (45): 


TES 1900 1% 
oí alluot) uot? 


Es claro que aparte de la primera derivada 
du/oít = du,/ 08 + du./0€, se requiere de la 
segunda derivada de las funciones u(5t) y 
u (Et). Éstas se obtienen de las ecuaciones (53) 
y (63), respectivamente: 


(65) 


uy MM - E 
f Aoc a 
dE 2 á ee E +) 


(66) 


O Once e 5 af 
qe 2/4 
6 E 
lar at 


=> ón, A esp[(a/96+1] ex E E " 08) 


=erfc 6, + 4) 


2 
op [- (4,/x)5 + ate E 


(siga) 


Lp [44/64 +(4,/0% 
+ AS) erfc Tb azi] 
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Figura 1. Las curvas de retención de humedad y de conductividad hidráulica de un suelo de Montecillo, México. 


Las características de Fujita y Parlange están definidas por las ecuaciones (11) y (12). 


Nótese que a partir de la segunda derivada 
espacial se puede calcular, considerando la 
ecuación (41), la derivada temporal: du/dt = 
(0%, 08 + 4.108). 


Aplicaciones 


La solución exacta es ilustrada en un suelo 
de Montecillo, México. De acuerdo con la 
clasificación textural americana, es un suelo 
franco arenoso. Este suelo fue caracterizado 
hidrodinámicamente con el método del drenaje 
interno por Cano (1990). 

En la figura 1 se muestran las características 
hidrodinámicas experimentales, así como las 
características de Fujita y Parlange, definidas 
por las ecuaciones (11) y (12), con parámetros: 
9, = 0.520 cm*/cm*, O, = 0.185 cm*/cm?, | y 
= 13.5 cm, | yw] =0 cm, K =25 cem/h, K_ =0 
cm/h, o. = 0.969, $ = 0.998. Se observa que las 
características teóricas representan de manera 
aceptable las características experimentales. 


En la figura 2, la solución es mostrada para 
una precipitación constante q, = 2 cm/h. En la 
parte (a) se muestra la evolución del perfil de 
humedad en una columna de suelo, con un 
contenido de humedad inicial constante 0, = 
0.250 cm*/cm? a todo lo largo de la misma. En 
la parte (b) se muestra la evolución del perfil 
de humedad, pero cuando la columna tiene 
dos contenidos de humedad, a saber 0, = 0.250 
em*/cm en 0 <z<z, y 6, = 0.500 cm*/cm!? en 
z,S<Z, donde z, =25 cm; la gráfica muestra que 
el agua de la zona más húmeda se distribuye 
hacia arriba por ascenso capilar y hacia abajo 
por efectos gravitacionales. 

Cuando el flujo adimensional de Darcy en 
la superficie de la columna tiene la expresión 
9. = 9, + yu, que resulta de la ecuación (34), la 
solución se muestra en las figuras 3 y 4. Con los 
valores de q,=2 cm/h y y=0.015, en la parte (a) 
de la figura 3 se tiene la evolución del perfil de 
humedad en una columna con un contenido 
de humedad inicial O, = 0.250 cm*/cm* y en 
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(a) 


(b) 


Figura 3. La solución exacta no lineal con término de sumidero, con parámetros q, =2 cm/h y y=0.015 y 
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Figura 2. La solución exacta no lineal, precipitación constante q, =2 cm/h. Humedad inicial: 
a) 0, = 0.250 cm*/cm', O < z; b) O, = 0.250 cm*/cm? en O <z <z, y 0, = 0.500 cm*/cm? en z, <z, donde z, = 25 cm. 
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condición inicial 0, = 0.250 cm? /cm?, O < z: a) evolución del contenido de humedad y b) evolución del término de sumidero. 
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Figura 4. La solución exacta no lineal con término de sumidero, con parámetros q,=2 cm/h y y=0.015, y 
condición inicial 0, = 0.250 cm? /cméen0<z< z, y 0, =0.500 cm?/ cm? en z ¡Sz: a) evolución del contenido 


de humedad y b) evolución del término de sumidero. 


la parte (b), la evolución correspondiente del 
término de sumidero. 

En la figura 4 se toma como condición 
inicial O, = 0.250 cm*/cm? en 0<z<z, y O, = 
0.500 cm*/cm? en z, < z, donde z, = 25 cm; la 
parte (a) muestra la evolución del contenido 
de humedad, mientras que la parte (b) muestra 
la evolución del término de sumidero. Se debe 
observar en este caso que el término desumidero 
toma tanto valores positivos como negativos 
debido a que es proporcional al gradiente del 
potencial de Kirchhoff y, en consecuencia, al 
gradiente de humedad, como se deduce de la 
ecuación (30). 


Conclusiones 


Se han planteado las ecuaciones diferenciales 
del flujo del agua en el suelo. Las trans- 


formaciones para resolver la ecuación de 
flujo han sido expuestas y desarrolladas, 
remarcando el significado físico de las 
mismas. De manera sistemática se han ex- 
puesto las características hidrodinámicas que 
han permitido la reducción de la ecuación 
Fokker-Planck de la infiltración a la ecuación 
de Burgers de la aerodinámica. Ellas han sido 
establecidas por Sanders et al. (1988) en la 
adaptación de la solución de la ecuación de 
Burgers a los medios porosos. 

La aplicabilidad de las transformaciones 
demanda características hidrodinámicas es- 
pecíficas, es decir, la difusividad de Fujita 
(1952) y la relación entre la conductividad y la 
difusividad propuesta por Parlange et al. (1982). 
Se demuestra que en el fenómeno del ascenso 
capilar, estas características proporcionan un 
volumen almacenado acotado en tiempos 
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grandes, a diferencia del volumen infinito que 
se obtiene con las características de Broadbridge 
y White (1988). 

Las mismas transformaciones han permi- 
tido la proposición de la variación posible 
del término de sumidero de la ecuación de 
Fokker-Planck. Las soluciones de la ecuación 
de Burgers han sugerido que el término de 
sumidero puede tener un comportamiento 
lineal en el potencial de Kirchhoff. 

La solución general aquí obtenida ha sido 
ilustrada en el suelo Montecillo, asimilando el 
término de sumidero a la extracción de agua 
por las raíces de las plantas. Por otra parte, 
la solución exacta construida puede ser útil 
en la validación de soluciones numéricas de 
la ecuación de transferencia del agua en los 
suelos. 
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Abstract 


FUENTES, C., CHÁVEZ, C., SAUCEDO, H. € ZAVALA, M. On an exact solution of 
the non-linear Fokker-Planck equation with sink term. Water Technology and Sciences, 
formerly Hydraulic engineering in Mexico (in Spanish). Vol. IL, No. 1, January-March, 
2011, pp. 117-132. 


The Fokker-Planck equation with a sink term is resolved exactly, using the Fujita diffusivity 
and Parlange relationship between conductivity and diffusivity. In order to obtain the solution, 
the Kirchhoff potential and the Fujita-Storm mobile coordinate are introduced. The differential 
equation takes the form of the Burgers equation, which is linear in the diffusive term. The 
convective coefficient of the latter is replaced by the Hopf-Cole transformation for the purpose 
of deriving the classical linear heat equation. During the transformation, the sink term is 
defined functionally, so that the end result is precisely the heat equation without sink term. The 
exact solution of the Hopf-Cole potential is deduced by using the classical Laplace transform 
for certain initial and boundary conditions of interest. The solution of the Fokker-Planck 
equation in the physical space is obtained through the inverse transformations. The solution 
includes as particular cases both the Sanders et al. and Broadbridge and White solutions. The 
exact solution can be used to validate numerical solutions of the Fokker-Planck equation and 
in studies on water extraction by plant roots. 


Keywords: Fujita-Parlange characteristics, Kirchhoff transformation, Hopf-Cole transfor- 
mation, Burgers equation. 
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